Modern meson— exchange potential and superfluid neutron star crust matter 
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In this work we study properties of neutron star crusts, where matter is expected to consist of 
nuclei surrounded by superfluid neutrons and a homogeneous background of relativistic electrons. 
The nuclei are disposed in a Coulomb lattice, and it is believed that the structure of the lattice 
influences considerably the specific heat of the neutronic matter inside the crust of a neutron star. 
Using a modern meson-exchange potential in the framework of a local-density approximation we 
calculate the neutronic specific heat accounting for various shapes of the Coulomb lattice, from 
spherical to non-spherical nuclear shapes. We find that a realistic nucleon-nucleon potential leads 
to a significant increase in the neutronic specific heat with respect to that obtained assuming a 
uniform neutron distribution. The increase is largest for the non-spherical phase of the crust. 
These results may have consequences for the thermal history of young neutron stars. 

PACS number(s): 97.60. Jd 2I.65.+f 74.25.Bt 



1 



The observation of thermal emission from the surface 
of a neutron star is a powerful tool by which one can 
obtain information about the state of matter inside the 
star. It has been shown that the time needed for a tem- 
perature drop in the core to affect the surface temper- 
ature should depend on the thickness of the crust and 
on its thermal properties, such as the total specific heat 
101, which is strongly influenced by the superfluid state 
of matter inside the crust. 

It has recently been proposed that the Coulomb-lattice 
structure of a neutron star crust may influence signifi- 
cantly the thermodynamical properties of the superfluid 
neutron gas The authors of Refs. ||-|^ have proposed 
that in the crust of a neutron star non-spherical nuclear 
shapes could be present at densities ranging from p — 
1.0 X 10^"* gcm~'^ to p = 1.5 X 10^'' gcm^'^, a density re- 
gion which represents about 20% of the whole crust. The 
saturation density of nuclear matter is po = 2.8 x 10^"* 
gcm"'^. These unusual shapes are supposed iQ,^ to be 
disposed in a Coulomb lattice embedded in an almost 
uniform background of relativistic electrons. According 
to the fact that the neutron drip point is supposed to 
occur at lower density {p ^ 4.3 x 10^^ gcm~^), and con- 
sidering the characteristics of the nuclear force in this 
density range, we expect these unusual nuclear shapes to 
be surrounded by a gas of superfluid neutrons. 

In the present paper we follow essentially Ref. 
however, we differ in using one of the modern meson- 
exchange potentials of the Bonn group Q in evaluating 
the neutron pairing energy gap. This potential is ex- 
pected to give, contrary to the effective one used in Ref. 
0, a more realistic estimate of the gaps in the region 
inside the nuclear cluster. 

In the following we treat the crust-lattice in the 
Wigner-Seitz approximation, dividing the Coulomb lat- 
tice into unit cells of appropriate shape (cylindrical, pla- 
nar and spherical), containing a nucleus surrounded by 
a gas of superfluid neutrons. Details of this system have 
been worked out within the framework of Thomas-Fermi 
calculations with different energy density functionals . 
The neutron and proton density profiles are obtained 
from a recent parametrization of Oyamatsu H 
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for r > Ri, where i = p,n represent protons and neu- 
trons respectively, with taken to be zero. Following 
Oyamatsu the spherical shape is expected for densities 
p < 0.35po (about 80% of the whole crust). The cylindri- 
cal region is supposed to appear for densities p between 
0.35po and 0.46po, while in the region OAGpo < p < 0.5po 
one presumes to have a slab-like form for the Coulomb 
lattice. The parameters Ri represent the finite bound- 



ary of the nucleus and ti determines the relative surface 
thickness. 

With these density profiles we already have a fixed 
proton fraction relevant for the calculation of the pairing 
gap. Our calculation of the pairing gap is a two-step 
process (for details, see Ref. j|]). First we solve self- 
consistently the Brueckner-Hartree-Fock (BHF) equa- 
tions for the single-particle energies, using a G-matrix 
defined through the Bethe-Brueckner-Goldstone equa- 
tion as 
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where V is the nucleon-nucleon potential, Q is the 
Pauli operator which prevents scattering into intermedi- 
ate states prohibited by the Pauli principle, Hq is the un- 
perturbed hamiltonian acting on the intermediate states 
and uj is the so-called starting energy, the unperturbed 
energy of the interacting states. Methods to solve this 
equation are reviewed in Ref. (|]. The single-particle 
energies for state ki (i encompasses all relevant quantum 
numbers like momentum, isospin projection, angular mo- 
mentum, spin etc.) in nuclear matter are assumed to have 
the simple quadratic form|^ 
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where m* is the effective mass. The terms m* and 6, 
the latter being an effective single-particle potential re- 
lated to the G-matrix, are obtained through the self- 
consistent BHF procedure. The so-called model-space 
BHF method for the single-particle spectrum has been 



used, see e.g. Refs. [^^|, with a cutoff fcM — 3.0 fm~ . 
This self-consistency scheme consists in choosing ade- 
quate initial values of the effective mass and S. The ob- 
tained G-matrix is in turn used to obtain new values 
for TO* and S. This procedure continues until these pa- 
rameters vary little. The BHF equations are solved for 
different proton fractions, using the formalism of Refs. 
1^,0. The nucleon-nucleon potential is defined by the 
parameters of the meson-exchange potential model of the 
Bonn group, version A in Table A. 2 of Ref. [Q. 

The next step is to evaluate the gap equation following 
the scheme proposed by Anderson and Morel [|l2| and 
applied to nuclear physics by Baldo et al. |jl^]. These 
authors introduced an effective interaction Vk^k' ■ This 
effective interaction sums up all two-particle excitations 
above the cutoff kM- It is defined according to 
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where the quasiparticle energy Ek is given by 
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Sf being the single-particle energy at the Fermi surface, 
Vk^k' is the free nucleon-nucleon potential in momentum 
space and is the pairing gap 

For notational economy, we have dropped the subscript 
i on the single-particle energies. 

In summary, first we obtain the self-consistent BHF 
single-particle spectrum Sk, thereafter we solve self- 
consistently Eqs. (j^) and (Q) in order to obtain the pair- 
ing gap A. This pairing gap is again calculated for var- 
ious proton fractions according to En. (1). For states 
above knj, the quasiparticle energy of (0) is approximated 
by Ek — {sk ~ ep), an approximation found to yield sat- 
isfactory results in neutron matter |q] . 



FIG. 2. Neutron effective mass ratio m*/m as function of 
the position r inside the Wigner-Seitz cell for the different 
densities and shapes reported in Fig. 1 



FIG. 1. Neutron pairing energy gap A as function of the 
position r inside the Wigner-Seitz cell for various nuclear 
shapes and densities. For the spherical phase, p/po = 0.058 
and p/po = 0.176 we have used solid and dashed lines, re- 
spectively. The cylindrical region (p/po = 0.354) is given by 
a dotted Une, while the slab region {p/po — 0.48) is shown 
with a dash-dotted line. 

The results for the neutron pairing gap, effective mass 
and local Fermi momentum as a function of the position 
in the Wigner-Seitz cell are displayed in Figs. 0, g and ||, 
respectively, for various density regions inside the crust. 



FIG. 3. Neutron Fermi momentum as function of the posi- 
tion r inside the Wigner-Seitz cell for various nuclear shapes 
and densities as reported in Fig. 1 

For heavy and medium heavy nuclei present in the neu- 
tron star crust, we expect mean radii of the order 6 ~ 8 
fm. We notice in Fig. ^ that the Fermi momentum for 
r < 6 fm is close to that of the saturation density of nu- 
clear matter, or the central density of ^"^Pb. At these 
densities, the gap energy in nuclear matter is generally 
small [Hjl^. This is also seen in Fig. |l|, where the gap 
is less than 0.5 MeV. The differences between the gap 
energies for the various shapes of the Coulomb lattice 
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can be retraced to the different effective masses in Fig. 

since the effective masses which enter the determina- 
tion of the pairing gap, differ. For values of r between 
6 and 8 fm, close to the Fermi surface, where the Fermi 
momentum of Fig. ^ changes rapidly, we see the largest 
variations in the effective mass and the pairing gap, the 
latter reaching a peak of 2.5 MeV for the spherical phase. 
For larger values of r, i.e. outside the nucleus, the Fermi 
momentum stabilizes, though there are large differences 
from shape to shape. For the spherical shape given by 
the lowest density p/po = 0.058 we get the lowest value of 
kp, whereas for the slab phase with p/ po — 0.48 (half nu- 
clear matter saturation density), we have kp « 1.3 fm^^. 
Since the neutron ^5*0 pairing gap reaches its maximum 
at kp between 0.8 ~ 1.0 fm~^ HJl^i we see from Fig. 
|l|, that the pairing gap is at its largest for the spherical 
phases and for radii larger than 8 fm. For large values of 
r, the pairing gaps are constant, since we have uniform 
neutron matter at a fixed kp. 

The qualitative features exhibited in Figs. |l|-^, are sim- 
ilar to the results of Broglia et al. |^. However, the 
pairing gaps obtained with the effective interaction of 
Ref. Ig], are larger than those obtained here. Broglia 
et al. obtain maximum pairing gaps of the order of 3.5 
MeV, whereas ours are of the order of 2.5 MeV. Our pair- 
ing g ap for uniform matter is close to that of Baldo et 
al. |13| , who also employ realistic nuclcon-nuclcon po- 
tentials. Moreover, in applications to finite nuclei the 
nucleon-nucleon force used here reproduces very well e.g. 
the experimental spacing between the 0^ ground state 
and the first excited 2+ state of the tin isotopes. The 
smaller energy gap in this work may in turn have impor- 
tant consequences for thermal properties of neutron stars. 
To see this, we evaluate the specific heat for a system of 
superfluid neutrons. Here we use the thcrmodynamical 
expression 
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temperature the system of superfluid neutrons exhibits a 
phase transition towards a normal Fermi liquid system. 
The contribution to the total neutronic specific heat in- 
side the Wigner-Seitz cell is written as 
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where Vnorm is the volume occupied by the normal sys- 
tem inside the cell. 

The total neutronic specific heat can be written for all 
temperatures as 
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In order to show the relevance of our results, we com- 
pute the total fermionic specific heat given by 



Ct — Cn + Ce 
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where Cg is the specific heat for relativistic electrons 

T 
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and ep IS the Fermi energy of the electrons ||lj]. The 
large anisotropy found in the local pairing energy gap 
(see Fig. ||), leads to a larger specific heat for the super- 
fluid neutron gas with respect to that obtained for the 
uniform neutron system. This is due to the weak super- 
fluid neutron component inside the nucleus, see Fig. 0. 
This effect yields a larger specific heat compared to that 
of the superfluid neutron matter outside the nucleus . 



where V is the volume of the system, T is the temper- 
ature and S is the total entropy. We obtain then the 
following form for the speciflc heat (see Ref. |^ for fur- 
ther details) 
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where E^p (r) is the local quasi-particle energy 

Eu, (r) ^ ^{e,{r)-epY + ^l^{r), (10) 

Vws is the volume of the Wigner-Seitz cell and Vsup is 
the volume occupied by the superfluid. With increasing 



FIG. 4. Ratio C^" /C^ as function of temperature T for 
spherical and non-spherical phases inside the crust. Notations 
as in Fig. 1 
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In Fig. ^ we compare the ratio CT{n.u.) / Ct{u.) be- 
tween the total fermionic specific heat evaluated account- 
ing for non-uniform nuclear shapes and that obtained 
considering uniform neutron matter only. We notice that 
the ratio increases dramatically moving from spherical 
to plate-like nuclei, and for T > 0.008 MeV, typical of 
the inner crust of a neutron star in the first 10^ yr af- 
ter formation, the non-spherical nuclear phases give spe- 
cific heats which are much larger than that for uniform 
neutron matter. This strong increase is due to the cor- 
responding increase in the nuclear volume observed for 
non-spherical phases [^6|. This increase overcomes the 
reduction in the anisotropy obtained in the pairing en- 
ergy gap moving from spheres to slabs. 

In summary, our results may have consequences for 
the thermal evolution of the star compared to models 
where the crust is described in terms of uniform neutron 
matter only. In fact properties like the heat diffusion 
time through the crust may be affected. Moreover, 
compared to the results reported by Brogha et al. [Q, 
where the Gogny interaction was employed to obtain the 
pairing gap, there is a further enhancement of the ratio 
CTin.u.)/CT{u.). 
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